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Abstract: Recent results concerning the topological properties of random geomet- 
rical sets have been successfully applied to the study of the morphology of clusters 
in percolation theory. This approach provides an alternative way of inspecting the 
critical behaviour of random systems in statistical mechanics. 

For the 2d g-states Potts model with q < 6, intensive and accurate numerics indi- 
cates that the average of the Euler characteristic (taken with respect to the Fortuin- 
Kasteleyn random cluster measure) is an order parameter of the phase transition. 
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1 Introduction 

Recently, new insights in the study of the critical properties of clusters in perco- 
lation theory have emerged based on ideas coming from mathematical morphology 
and integral geometry [Ha, ^, 3cE|. These mathematical theories provide a 



set of geometrical and topological measures allowing to quantify the morphological 
properties of random systems. In particular these tools have been applied to the 
study of random cluster configurations in percolation theory and statistical physics 
PW| , 0, |, [Wall , [Ml . 



One of these measures is the Euler-Poincare characteristic x which is a well 
known descriptor of the topological features of geometric patterns. It belongs to 
the finite set of Minkowski functionals whose origin lies in the mathematical study 
of convex bodies and integral geometry (see P4 B |S c!t|] ) . These measures, as we 
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shall explain below, share the following remarkable property: any homogeneous, 
additive, isometry-invariant and conditionally continuous functional on a compact 
subset of the Euclidean space M. d can be expressed as a linear combination of the 
Minkowski functionals. This is the well known Hadwiger's theorem |[Ha|| of integral 
geometry which has a wide scope of applications in mathematical physics due its 
rather general settings. 

The use of these measures in image analysis ||Gra|| , problems of shape recognition 
3e], determination of the large scale structures of the universe |[Wa2|| , modelling of 



porous media [ AKPM |, microemulsions | HK | and fractal analysis [|M2 | has been a 
topic of growing interest recently. 

The application of these description tools for the study of random systems in 
statistical mechanics has already provided interesting results. In |[MW| , Q, the 
computation of the Euler-Poincare characteristic \ f° r a system of penetrable disks 
in several models of continuum percolation has led to conjecture new bounds for the 
critical value of the continuum percolation density. In [ MW |, an exact calculation 
of x has shown that a close relation exists between the zero of the Euler-Poincare 
characteristic and the critical threshold for continuum percolation in dimensions 2 
and 3. 

Of similar interest in the same domain, let us recall that for the problem of bond 
percolation on regular lattices, Sykes and Essam |[SyEs|| were able to show, using 
standard planar duality arguments, that for the case of self-dual matching lattices 
(e.g. Z 2 ), the mean value of the Euler-Poincare characteristic changes sign at the 
critical point (this even led them to announce a proof for the existence of the critical 



probability of bond percolation on Z 2 ), see also |Gri 



More recently, Wagner |[Wal| | was able to compute the Euler-Poincare charac- 
teristic on the set of all plane regular mosaics (the 11 archimedean lattices) as a 
function of the site occupancy probability p G [0, 1]] and showed that a close con- 
nection exists between the threshold for site percolation on these lattices and the 
point where the Euler-Poincare characteristic (expressed as a function of p) changes 
sign. 

The aim of this article is to further investigate the role played by this mor- 
phological indicator in statistical physics and to present new results concerning its 
behaviour in the case of the 2-dimensional Potts model. Namely we present here 
clear evidence, based on Monte Carlo simulations, that for the 2d Potts model, the 
Euler-Poincare characteristic is an order parameter of the phase transition. Namely 
we find that for q — 2, 4 it changes sign continuously at the transition point while, 
for g = 5,6 it has a first order transition at the critical point. As far as we know, 
this is the first example of a discontinuous behaviour of this parameter in a physical 
model. 

The paper is organized as follows. In Section 2 we introduce a minimal back- 
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ground concerning Minkowski functionals and the necessary definitions for our model, 
the numerical results are presented in Section 3 followed by some comments and dis- 
cussion in Section 4. 



2 The model 



We first briefly summarize the basic facts from integral geometry and give the defi- 
nition of the Minkowski functionals including the Euler-Poincare characteristic, see 
Ha| , ^, |Sch|] for more complete expositions. We will then show how to compute the 



Euler-Poincare characteristic in the case of a random configuration of sites and bonds 
produced by a Fortuin-Kasteleyn (FK) transformation of the partition function of 
the Potts model [IFK . 



Minkowski Functionals 

In the most abstract setting, the Minkowski functionals are defined as integrals 
of curvatures in the framework of the differential geometry of smooth surfaces in 
compact sub-domains of the Euclidean space 

Another way to introduce these measures is as follows. First define the Euler- 
Poincare characteristic x as an additive functional on subsets of H d such that for 
A, B C H d (see example in Fig. 1) 



X {A UB) = X {A) + X (B) - x{A n B) 



with 



1, if A ^ 0, A convex 

o, A = m. 



(2.1) 




Figure 1: Examples of calculation of the Euler-Poincare characteristic in dimension 
d = 2 for some combinations of convex subsets of H 2 . 
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The set of Minkowski functionals is then defined by (see e.g. ||Gra| , [M\V| ) 

W a (A) = j 'x(AUE a )dti{E a ), a = 0,...,d-l (2.2) 

W d (A) = co dX (A): co d = n d / 2 /T(l + d/2) (2.3) 

where E a is an a-dimensional plane in M d , dfj,(E a ) its density normalized such that, 

and ujd is the 



d— a 



for the d-dimensional ball B d (r) with radius r, W a (B d (r)) = tu d r 
volume of the unit ball in IR d . 

Obviously, additivity of the Minkowski functionals is inherited from (|2.1|) , fur- 
thermore they are conveniently normalized through 



M a (A) 



LO d LO a 



-W a (A) 



The computation of these normalized functionals in dimensions 1, 2, 3 in term of the 
usual geometric measures (length, area, volume,...) is given in Table 1. 



d 
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M 3 
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Xi/2 
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H/2-K 2 


3x 3 /4tt 



Table 1: Values of the normalized Minkowski functionals for <i-dimensional subsets 
of !R d , d — 1, 2, 3 in terms of geometric measures, L: length, S: area, V: volume, C: 
circumference, H: integral mean curvature, \ d . Euler-Poincare characteristic. 



Hadwiger's completeness theorem asserts that, under not too restrictive and fur- 
thermore physically reasonable assumptions, namely: additivity, motion invariance 
(under translations and rotations) and conditional continuity (which states that any 
convex body can be smoothly approximated by convex polyhedra), any functional 
M. decomposes as a linear combination of the finite set of Minkowski functionals 

d 

M{A) = J2 c » M v( A ) 

where the c v are real coefficients. 

This theorem has very important practical consequences, we refer to |[Ha| , |^, |Gra 
we shall only mention the 'principal kinematic formula (see P, Mil) 



I M a (AUgB)dg = ^( C *)M a ^(B)Mp{A), a = 0,...,d 
Jg p =0 \P/ 
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where integration is over the group of motions (i.e. rotations and translations, 
g = (r, ©)). This formula is very useful to calculate mean values of Minkowski 
functionals for random distributions of objects. For instance, it can be applied 
to the computation of the excluded volume of convex bodies (leading, in case of 
spherical objects, to Steiner's formula), which has been used to estimate the critical 
thresholds in continuum percolation theory |[Ba| , |FS . 



In the following, we will be mainly interested in the evaluation of the Euler- 
Poincare characteristic for random bond configurations on regular lattices A C 7L d 
which, as we shall see, is equivalent to compute Euler formula for planar graphs 
made of S sites (or nodes), B bonds (or links) and V plaquettes (or faces) ||Agj , [Be] , 
i.e. 

X =S-B+V 

It is indeed one of the topics of algebraic topology to show that the above general 
definitions of Minkowski functionals extend to cell complexes (see Appendix) to 
which random bond configurations on regular lattices belong. 

Because of the growing interest in the use of these morphological measures, we 
thought that it was important to recall briefly their definitions and main properties. 



Potts Model 

The partition function for the g-states Potts model on A C Z d at inverse temperature 
(3 reads 

^ p r s ( A ) = E ex p{' 9 E <Wi)} ( 2 - 4 ) 

where the first sum runs over all configurations er C {1, . . . , g}' A ', the second one is 
over each nearest neighbour pair of Potts spins on A and 5 is the Kronecker symbol. 
We remind that, whenever q is large enough, in any dimension d > 2, this model 
exhibits a unique (inverse) temperature f3 c where the mean energy is discontinuous 
(see IKl [LMP4 [KLMRp . In dimension d = 2 for q > 5 this is an exact result |[Ba:x 



and it is expected to be true, in d = 3, for q > 3 [|Wu 



After performing the (FK) transformation [ FK , the partition function ( |2.4| ) leads 



to the following random cluster representation 

ZJJ(A) = - l)»Hx) q mx) (2.5) 



x 



Here the summation is over all graphs X which can be drawn inside the domain A, 
N l (X) is the number of bonds of the configuration X and N\(X) is the number 
of connected components of X (including isolated sites). We will call A^°(X) the 
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number of sites which are endpoints of a bond in X and N 2 (X) the number of 
plaquettes of the configuration X, i.e. the set of cells in A having 4 occupied bonds 
on its boundary (see Fig. 2). 
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Figure 2: A bond configuration X C B\ with 25 sites, 19 bonds and 2 plaquettes. 



In this framework, the Euler-Poincare formula ||Ag|| leads to the following defini- 
tion of the Euler-Poincare characteristic 

X (X) = N°(X)-N 1 (X) + N 2 (X) (2.6) 

This expression will allow us to compute the mean value of the Euler-Poincare 
characteristic with respect to the FK measure. 

In algebraic topology, one is interested in classifying geometric objects (or spaces) 
up to some transformations acting on these spaces (homeomorphisms) by means of 
simple invariant quantities. It turns out that this program can be achieved and such 
invariants of a space are the homology groups of these spaces. In this framework, the 
Euler-Poincare characteristic (as the other Minkowski functionals) can be expressed 
in terms of these invariants and lead to a definition for x equivalent to the one given 
by (|2.6|) (see Appendix). 
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3 Numerical Results 

We have performed Monte Carlo simulations of the 2-dimensional g-state Potts 
model for q ranging from 2 to 6. We have always simulated the models near the 
critical temperature, whose value can be exactly determined through the well known 
formula: (3 c (q) = J/kT c (q) = log(l + y/q). In order to extract a value of the Euler 
characteristic \ as close as possible to the value at the infinite volume limit, we 
have taken rather large lattices: for the three models with a continuous transition 
(q = 2,3,4) we arrived at lattice sizes up to 2000 2 . The algorithm we used is 
the Wolff cluster update |[Wo |. The identification of FK cluster configurations has 



been performed via the Hoshen-Kopelman algorithm [[Ho|| ; we always considered 
free boundary conditions for the cluster labeling. 

The Wolff algorithm is the less efficient the bigger the number q of states. Par- 
ticularly dramatic is what happens when one passes from the 2-state (Ising) to the 
3-state model: in the former case, on the 2000 2 lattice it is enough to perform few 
updates (5 — 10) to get uncorrelated configurations for the cluster variables, in the 
latter one needs about 1000 updates! Because of that, the simulations for q = 3,4 
on the 2000 2 lattice were very slow, and the relative data could not reach a high 
statistics. Nevertheless, as we will see, we can get useful indications out of them. 

If x changes sign at the threshold, on a finite lattice the values measured at each 
iteration would be distributed around zero, provided the lattice is large enough. 
Therefore one would see both positive and negative values. For this reason, it is 
helpful to look at the distribution of x- I n the following we present separately the 
results for q = 2, 3, 4 and q — 5,6. 



3.1 Results for the models with a continuous phase transi- 
tion 

The first case we consider here is the Ising model. Fig. 3 shows the x distribution for 
three different lattice sizes: 500 2 , 1000 2 and 2000 2 . In each case we have taken 20000 
measurements. The peak of the distribution shifts towards x = the larger the 
lattice. The average values of x are: x(500 2 ) = 0.00101(2), x(1000 2 ) = 0.00053(2), 
x(2000 2 ) = 0.00024(1). We notice that the averages are quite small and decrease 
sensibly if we go to larger sizes, reducing themselves to about the half when we 
pass from a lattice to the next one. This approximate linear scaling of % with the 
lattice side suggests that the Euler characteristic at the infinite volume limit indeed 
vanishes. 

Let us now examine the case q = 3. In Fig. 4 we again plot the x distribution 
for the same three lattice sizes we have considered for the Ising model. 
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2500 




Euler characteristic X 

Figure 3: Distribution of x fc> r the FK cluster configurations of the 2d Ising model 
at the critical point. 

900 i 1 1 1 -r 1 1 1 1 1 




Euler characteristic X 

Figure 4: Distribution of x f° r the FK configurations of the 2d, 3-state Potts model 
at the critical point. 
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Since we collected a different number of measurements for the different lattices, 
for a real comparison of the distributions we needed to renormalize the total number 
of measurements on each lattice to the same value: we decided to renormalize all data 
sets to the number of measurements on the 500 2 lattice (10000). The distributions 
are broader than the Ising ones but they appear almost exactly centered at x — 0- 
The average values are in fact much smaller than before: x(500 2 ) = 0.00024(4), 
x (1000 2 ) = 0.00012(3), x(2000 2 ) is zero within errors. We then deduce that also for 
g = 3 x = at the critical point. 

To complete our analysis we studied the case q — 4. In Fig. 5 we present a 
comparison of the x distributions for two lattice sizes, 800 2 and 2000 2 . 




Figure 5: Distribution of x fc> r the FK cluster configurations of the 2d, 4-state Potts 
model at the critical point. 



There is a clear shift of the center of the distribution towards zero when one goes 
from the smaller to the larger lattice. The average values of the Euler characteristic 
in the two cases are x(800 2 ) = -0.00145(8) and x(2000 2 ) = -0.0006(2). Also here 
there is no apparent convergence to some value, even if the lattices are rather large: 
| < x > I reduces itself to less than its half by changing the lattice size. From all this 
we also deduce that the Euler characteristic of the FK clusters of the 2-dimensional 
4-state Potts model vanishes at criticality. 
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3.2 Results for the models with a first order phase transition 

For q > 4 the 2d g-state Potts model undergoes a first order phase transition, i.e. the 
thermal variables vary discontinuously at the critical threshold. The magnetization, 
for instance, makes a jump, varying from zero to a non-zero value. Because of that, 
we expect that the cluster configurations change abruptly at the critical point, and 
that the cluster variables exhibit as well discontinuities. In particular, the Euler 
characteristic may jump from a value to another. 

We analyze here the 5- and 6-state Potts models. In both cases we have per- 
formed simulations on a 300 2 lattice. In Figs. 6 and 7 we compare the distribution 
histograms of the magnetization M and the Euler characteristic x a ^ three different 
temperatures: above, near and below T c . We define the magnetization by taking 
the excess of sites in the majority spin state with respect to the value 1/q in the 
paramagnetic phase, when all spin states are equally distributed. Therefore we al- 
ways measure M > 0. Looking at the magnetization histograms one clearly sees the 
spontaneous symmetry breaking by reducing the temperature. The characteristic 
double peak structure of M around T c indicates that the transition is first order, as 
it is known. The corresponding histograms of the Euler characteristic show a per- 
fectly analogous pattern. As we expected, the two coexisting FK-phases at T c are 
characterized by two different values of the Euler characteristic (double peak in the 
figures). The ordered phase is made of bonds between Potts spins of the same color 
at T c . We then deduce that x, like M, varies discontinuously at the threshold. The 
result is valid for the 5-state and the 6-state Potts model, so it is likely to be valid 
also for q > 6, when the discontinuity of M at the threshold is sharper. Looking 
at both figures we remark that the centers of the peaks of x look approximately 
symmetric with respect to zero. If this symmetry exists, it would be an interesting 
feature, and at the moment we have no arguments to justify it. In order to determine 
with some accuracy the values of x i n the two coexisting phases we would need to 
increase considerably the size of the lattice, but the required computer time would 
increase dramatically for the reasons we explained at the beginning of this section. 
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Figure 6: Distribution histograms of the magnetization M and the Euler charac- 
teristic x f° r the 2- dimensional 5-state Potts model at three different temperatures. 
The lattice size is 300 2 . The behaviour of x is driven by M. 
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Figure 7: Distribution histograms of the magnetization M and the Euler charac- 
teristic x f° r the 2- dimensional 6-state Potts model at three different temperatures. 
The lattice size is 300 2 . The behaviour of x is driven by M. 
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4 Conclusions 

This work clearly indicates and confirms that the Euler-Poincare characteristic is 
indeed an important indicator of a phase transition, playing the role of an order 
parameter for the 2-d Potts model, to the extend of the cases (q = 2, 6) studied 
here. For this model, it reveals that the topology of cluster configurations has a 
deep meaning concerning criticality. 

The fact that x changes sign at T c can be understood in the following way. Let 
us consider the 2D Ising model: from T c up to T = oo the system is in its disordered 
phase and the only excitations one can get in the FK-bond representation are made 
of isolated bonds (the probability to see any plaquette vanishes exponentially). Ap- 
plying ( |6.1| ) (see Appendix), one sees that \ behaves like /3°(X) times a term of the 
order of the volume of the system. However, from T c down to T = 0, the system 
is in its ordered phase and the corresponding FK-configuration is (with high prob- 
ability) made of 0(1) connected bond components. Missing bonds constitute the 
excitations and their number scales with the volume of the system so, using again 
( p.l| ), one gets that \ behaves like —/3 l (X) times a term of the order of the volume 
of the system. This explains in the case of the Ising model the change of sign of x 
at T c . 

Other spin systems have to be investigated in order to see whether this property 
of the Euler-Poincare characteristic is shared by models with continuous symmetries 
such as the X-y-model or the Widom-Rowlinson model. 

Another important question concerns the critical behaviour in gauge models. 
For example a similar study could provide some insights concerning the deconfining 
transition in SU(N) gauge theory. Indeed, some works fSaj tend to indicate that this 



transition could be probed by percolation of some physical clusters related to color 
fields in lattice QCD. These models have been thoroughly investigated in the past 
and the tools coming from algebraic topology have been of primary importance to 
uncover profound duality results concerning their phase structure ||We| , |Wa3| . 



The interest of the last remark resides in the fact that so far, no suitable order 
parameter is known that can describe the deconfining transition in SU(N) gauge 
theory in the case of finite quark masses and it would be interesting to figure out 
how does the Euler-Poincare characteristic behave in gauge models. 
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6 Appendix 

The Euler-Poincare characteristic has another equivalent formulation showing closer 
connections with the topological properties of configurations. FK bond configura- 
tions X, together with a suitably defined orientation, can be endowed with a struc- 
ture of 2- dimensional (closed) cell-complexes (see [TA], |Wa3| , [KLMR| ). Sites, bonds 



and plaquettes are then called, respectively, 0-cells, 1-cells and 2-cells. A cell com- 
plex is then defined as a finite collection of g-cells, q = 0, 1, . . . 

A theorem in Algebraic topology states that (see e.g. [|Agi| ): the only topo- 



logically invariant functions on cell complexes are those which are functions of the 
Euler-Poincare characteristic and the dimension of the complexes. 

Formal linear combinations of oriented g-cells are called g-chains and form an 
additive group. Of primary importance is the definition of a boundary map d q 
acting on a g-cell (more generally on g-chains) giving rise to a (q — l)-cell, q > 1 
(the boundary map of a 0-cell being defined as a 0-cell). For example the boundary 
map acting on an oriented plaquette (2-cell) leads to the set of all oriented bonds 
(1-cells) defining its boundary. 

Let a g-cell. If d q c^ q ' = then c- q > is called a q-cycle and if S q+1 ^ is a 
(q + l)-cell such that d q+ i c^ +1 ^ = c^ q > then c^ +1 ^ is called a q-boundary. The set 
Z q of g-cycles and B q of g-boundaries form subgroups of the group of g-chains. The 
quotient group H q = Z q /B q is called the g-th homology group and it is the topic of 
algebraic topology to show that the rank of these groups (called the Betti numbers) 
is invariant under homotopy transformations. 

This means, loosely speaking, that there is a way to identify geometric struc- 
tures through continuous deformations without introducing or suppressing "holes' 
in them; for example a sphere in 1R 3 is not homotopy equivalent to a point. In 
our context of clusters of bond configurations, a connected component (collection 
of plaquettes and bonds) can be shrinked (and so it is homotopy equivalent) to a 
point. In that case it turns out that the rank of the 0-th homology group (the 0-th 
Betti number) has a very simple meaning, it just gives the number of connected 
components of the configuration. 

If j3 l denotes the i-th Betti number, then for two dimensional cell complexes like 
the one we consider in this paper, one has 

X {X)=(5\X)-(3\X) (6.1) 

Here f3 1 {X) (the rank of the first homology group) is the number of holes in the 
configuration, i.e. cycles which do not surround a connected collection of plaquettes, 
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each of them having 4 bonds on their boundary (see Fig. 8). One can verify that 
Eq. ( |6.1| ) applies to all examples given in Fig. 1. 

In dimension d = 3, /3 2 (X) would give the number of cavities inside the cell 
complex. 



a cycle that 

"bound" 
2 plaquettes 



r 



a "hole" : 
cycle that does not 
bound" anything 



Figure 8: Cycle and hole in FK-configuration. 



Equation (|6.1| ) shows that the Euler-Poincare characteristic of a bond config- 
uration depends only on homotopy equivalence classes of objects and as such, it 
is an intrinsic property which provides a distinguished signature of its topological 
features. The fact that it allows to probe the critical behaviour of physical systems 
makes possible the interpretation of phase transition on purely topological grounds. 

This formalism of cell complexes has been shown to be very useful in order to 
define general duality arguments in several models of statistical mechanics and it is 



the basis of the algebraic approach of the study of phase transitions [Wa3 
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